The Spontaneous Generation of Magnetic Moment of Beryllium Dimer on Graphene 
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Graphene has various potential applications in electronics and scientists are seeking to introduce 
magnetic properties on graphene. Here we report our theoretical findings that a local magnetic 
moment as large as 1 [ib can be generated from the adsorption of the diamagnetic beryllium dimer 
on perfect diamagnetic graphene. Unexpectedly the adsorption system is remarkably stable as 
evidenced by the binding energy of 1.04 eV, indicating that the formation of beryllium dimer is 
energetically much favorable. Analyses reveal that the magnetic moment and stability is originated 
from the spontaneous transfer of one electron from the anti-bonding orbital of the dimer to graphene. 

PACS numbers: 73.20.-r, 73.22.Pr, 75.75.Lf 
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Graphene has become a promising material for future 
electronics. The special electronic and magnetic proper- 
ties of graphene arc related to its low-dimensional carbon 
structures and have aroused enormous interest. [J, 0, Q 
Since the pure graphene is diamagnetic, introducing 
magnetic properties to graphene-related systems has at- 
tracted theoretical concern. 

H I, H By doping mag- 
netic transition-metal elements on perfect graphene or 
graphene vacancies one can introduce magnetic moment 
into graphene system. @, d, Adsorbing small covalent 
molecules on perfect graphene can also change the total 
magnetic moment by charge transfer between adsorbates 
and graphene. [13, El E3] It is notable that these cova- 
lent molecules stay away from graphene layer by 3 A at 
least and therefore the binding energies are small. Typi- 
cally, the transferred charges are less than 0.1 e~. Thus, 
the induced change of magnetic moment is expectedly 
small. 0, EH it is highly desired to find an adsorbate 
that can adsorb on graphene firmly and induce a large 
charge transfer or a remarkable magnetic moment. With 
density functional theory calculations, we found a large 
local magnetic moment can be produced spontaneously 
by adsorbing diamagnetic beryllium dimer on perfect dia- 
magnetic graphene. Here we report this interesting find- 
ing. 

All the calculations were performed with the plane 
wave based VASP code [H, QtSJlil using the projector- 
augmented wave methodjl*7l, Il8| and Perdew-Burke- 
Ernzernhof generalized gradient approximation [lEj. The 
supercell is composed of 4 x 4 unit cells with a distance 
of 15 A between the adjacent graphene layers. For the 
integration over the Brillouin zone, we combined 5x5x1 
Monkhorst-Pack grids [2^1 with a generalized Methfessel- 
Paxton smearing technique 2l| to optimize adsorption ge- 
ometries. Charge transfers were calculated based on the 
Bader charge analysis. 22, [23| For accurate charge and 
energy calculations a more fine 11 x 11 x 1 Monkhorst- 



Pack grids with tetrahedron method with Blochl correc- 
tions were used. An energy cutoff of 600 eV was used 
throughout the calculations. 

Beryllium dimer is the simplest metal dimer in ad- 
dition to lithium dimer. It has eight electrons totally 
and the electronic configuration of the ground state is 
(ler) 2 (lcr*) 2 (2cr) 2 (2cr*) 2 . The ground state of Be 2 is sin- 
glet since its electronic configuration is a type of closed 
shell. In old-fashioned text book, Be2 molecule cannot 
exist in reality because the full-filled anti-bonding 2a* 
orbital counteracts the full-filled bonding 2a orbital en- 
ergetically. Since 1980s, theoretical studies gradually in- 
dicate the existence of the stable Be2 with a short bond 
length by mixing excited determinants. 0, 25 1 Recent 
experimental work determined potential energy curve of 
Be2 and found it looks like a normal covalent molecule of 
weak bond strength with a bond distance of 2.45 A.[26j 
Our calculated bond length of beryllium dimer is 2.43 A, 
in nice agreement with the experimental value. 

Fig. [lja) shows the band structure of pure graphene. 
As comparison, the relative energy level of isolated Be 
2s, 2p atomic orbitals and Be2 2a, 2a* , l7r molecular or- 
bitals are also presented in Fig. QJb). The full-filled 2a* 
orbital is 0.33 eV lower than the Dirac point of graphene 
while the empty l7r orbital is 1.30 eV higher than the 
Dirac point. Such energy level ordering should prevent 
charge transferring from Be2 to graphene or vice versa. 
In fact, in the physisorption processes of diamagnetic co- 
valent molecules, such as H2O, NH3 and CO, the highest 
occupied molecular orbital (HOMO) and lowest occupied 
molecular orbital (LUMO) of the adsorption system are 
separated by the Fermi level so the total charge transfer 
is small. [lit 

Now we study beryllium dimer adsorbed on graphene. 
Our calculation revealed that though one single beryl- 
lium atom can hardly adsorb on graphene, the beryl- 
lium dimer does adsorb on graphene. We have investi- 
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FIG. 1: (a) The band structure of graphene. (b) The en- 
ergy levels of Be 2s, 2p atomic orbitals and Be2 2a, 2a* , In 
molecular orbitals. 



gated several adsorption configurations. The most sta- 
ble adsorption geometry is shown in Fig. O The Be- 
Be bond length on graphene is shortened by 0.32 A 
compared with that in gas phase. The perpendicular 
distance to graphene layer is 1.49 A (Fig. [5]) which 
is significantly shorter than those for other molecules 
doped on graphene. fiol El [H, H3| The binding energy 
(BE) of the dimer is unexpectedly as high as 1.04 eV. 
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large binding energy, it is also amazing that when the sin- 
glet diamagnetic beryllium dimer adsorbs on diamagnetic 
graphene, the Be2-graphene adsorption system acquires 
magnetic moment of 1 \ib per supercell. The emergence 
of magnetic moment usually indicates selective charge 
transfer of spin channel, which is contrary to our pre- 
vious prediction based on the energy spectra of isolated 
beryllium dimer and graphene. 




FIG. 2: The stable adsorption of beryllium dimer on the 
(4 x 4) cell of graphene. The dimer is perpendicular to the 
graphene and adsorbs on the hollow site. Be a atom is below 
Bet atom. Labeling of atoms: C, gray spheres; Be, yellow 
spheres. 

The band structure of Be2-graphene system is illus- 
trated in Fig. [31 The upshift of the Fermi level relative 
to the Dirac point results from the charge transfer from 



Be2 molecule to graphene. One can also find the split- 
ting of 2a and 2cr* spin orbitals of Be2 (see Fig. [3]). 
The 2cr spin orbitals are still lower than the Dirac point 
compared with the one in gas phase, but the 2a* spin or- 
bitals are pushed higher than the Dirac point. Fig. Ela) 
shows the spin-up component of 2a* ortital is below the 
Fermi level, while the spin-down component of 2a* orti- 
tal in Fig. [2[b) is higher than the Fermi level. Especially, 
the spin-down component of 2a* orbital is much higher 
than the Dirac point. Thus it is the electron of spin- 
down component on 2a* orbital which is transferred to 
the conduction band of graphene spontaneously. In Fig. 
IHa) we plot the total density of states (DOS) of two 
beryllium atoms and the total DOS of six carbon atoms 
around the hollow site. It is not difficult to verify that 
the split of 2a spin ortitals is about 0.35 eV, the 2a* or- 
bital of spin-up component is lower than the Fermi level 
by 0.27 eV and the 2a* orbital of spin-down component 
is higher than the Fermi level by 1.77 eV. Using Bader 
charge analysis we found three valance electrons stay on 
Be2 and the total magnetic moment of 1 \xb localizes on 
Be2 in the Be2-graphene system which means one elec- 
tron of spin-down component is transferred to graphene 
totally. 




FIG. 3: The band structure of Be2-graphene system, (a) 
The spin-up component, (b) The spin-down component. The 
2<T and 2a* orbitals of Be2 are marked in green and yellow, 
respectively. 



It is instructive to compare the adsorption of paramag- 
netic NO2 and diamagnetic Be2 on graphene. The mech- 
anisms of the charge transfer are different for the two 
adsorbates. In the case of the adsorption of the param- 
agnetic NO2 on graphene, it was argued that since the 
empty partially occupied molecular orbital (POMO) of 
adsorbed NO2 is below the Dirac point in the adsorbed 
configuration, the char ge t ransfer can reach one electron 
in the dilute limit. 13, 12, [28| In the actual finite theo- 
retical calculation, the empty POMO is crossed by the 
Fermi level so just fractional charge is transferred to the 
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FIG. 4: Comparison of local DOS and COOP, (a) Local DOS 
of two beryllium atoms and six carbon atoms around the ad- 
sorption sites, (b) COOP in red is of two beryllium atoms. 
The black one is the total COOP between Be a (see the cap- 
tion of Fig. [2} and six carbon atoms around the adsorption 
site. 



POMO from graphene and the total magnetic moment 
decreases. 13, U| In the Be2-graphene system, the Fermi 
level locates between the two distinctly split 2a* spin or- 
bitals of Be2 in the adsorbed configuration. Hence, one 
electron is transferred from the spin orbital above the 
Fermi level to graphene and a magnetic moment of 1 [Lb 
is produced spontaneously. 

Large binding energy of doping adsorbates on graphene 
usually involves the strong interaction between adsorbate 
atoms and carbon atoms of graphene, such as the adsorp- 
tion of transition-metal elements. Q To check whether 
this is the case in the Be2-graphene system, we also plot 
the total crystal orbital occupation populations (COOP) 
of adsorbed Be2 and the one between Be2 and six carbon 
atoms around the adsorption site in Fig. 0Jb). COOP 
is a more illustrative scheme to reflect the nature and 
strength of bonding or anti-bonding interaction between 
two bonding atoms than DOS.[29| The COOP plots are 
characterized with 2a and 2a* of Be2 orbitals. The 2a* 
spin orbitals are split by 2.04 eV which coincides with 
the DOS analysis. From Fig. 0Jb) we can infer that the 
interaction of Be2 2a* orbitals with graphene carbons 
is anti-bonding and negligible. The adsorption of beryl- 
lium dimer on graphene is a kind of physisorption essen- 
tially. But why the binding energy calculated is much 
larger than that of other small covalent molecule dop- 
ing on graphene where the binding energies is an order 
of meV? This is because the electron transferred from 
Be2 to graphene is from the anti-bonding 2a* orbital. 
The removal of one electron from the anti-bonding or- 
bital stabilizes Be2 and yields large binding energy. The 
shortening of the Be-Be bond length by 0.32 A provides 
a strong support for this point of view. 



By far, the study presented here only considers the 
ferromagnetic configuration of the Be2-graphene system. 
In order to check the magnetic state of the adsorption 
structure, we doubled the geometry in one direction and 
set the initial magnetic moments of beryllium to be an- 
tiferromagnetic. Calculation results showed the energies 
of both configurations are nearly equal, which means the 
size of supercell in our calculations is large enough to 
omit the magnetic couplings between Be2 of each super- 
cell. In addition, the spin-unpolarizcd configuration was 
also considered and found unfavorable energetically. 

In summary, we found beryllium dimer is favorable 
to form on graphene. The physisorption of diamagnetic 
beryllium dimer on pure graphene induces a charge trans- 
fer of one electron and generates a magnetic moment of 
1 fiB- The magnitude of charge transfer and magnetic 
moment generated by doping beryllium dimer are larger 
than by doping those reported molecules on graphene. 
Our study demonstrates that even without transition- 
metal adatoms or defective graphene a large magnetic 
moment can be generated spontaneously from two dia- 
magnetic matters. This opens up the possibility of intro- 
ducing stable and observable magnetic properties, such 
as magnetic domain and magnetic order, in graphene by 
molecular doping directly. 
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